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1. Introduction 

Over years significant efforts have been made to use lattice teciiniques to investigate the struc- 
ture of the nucleon. Of particular interest are the Parton Distribution Functions (PDFs) and form 
factors. The latter encode information about charge distribution and magnetization while the PDFs 
tell us about the distribution of momentum and spin. While some of the related observables can be 
determined with good accuracy by experiments (e.g. the nucleon's axial charge ^a) other quantities 
are difficult to access (like the tensor charge gx)- 

A precise determination of moments of nucleon PDFs and form factors on the lattice turned 
out to be rather challenging. It continues to be difficult to reach sufficient control on all system- 
atic errors such as finite size effects, lattice artefacts and the influence of the chiral extrapolation. 
Simulations are performed in volumes of a size where some quantities exhibit significant finite size 
effects. The available lattice data of the quantities of interest show no significant discretization ef- 
fects. But current simulations only probe a small window of lattice spacings thus providing us with 
limited control on the continuum extrapolation. From chiral effective theories (ChEFT) there are 
indications that the quark mass dependence close to the physical pion mass is very strong. There- 
fore extending lattice simulations into the region where rtiji < mps < 300 MeV has become a major 
goal for recent calculations. 

2. Simulation details 

For our simulations we use Wilson glue and N{ = 2 degenerate flavours of Clover fermions, 
where the improvement coefficient csw has been determined non-perturbatively. Most of our con- 
figurations have been generated using the BQCD implementation of the HMC algorithm Var- 
ious algorithmic improvements have been applied which accelerate this algorithm, such as the 
Hasenbusch preconditioning and the use of multiple time-scales, or reducing the time spent for 
matrix inversion. For instance, chronological guess and the Schwarz Alternating Procedure (SAP) 
are used to start the inversion and to precondition the fermion matrix 

These algorithmic improvements plus recent increase in computing resources enabled simula- 
tions in the region of small quark masses, i.e. in a region where the pseudo-scalar mass is smaller 
than 300 MeV. Fig. |l^ shows the parameter region of our simulations. When approaching physical 
quark masses larger lattices are needed to stay in the region mpsL > 3 where finite size effects are 
expected to be sufficiently small (see Fig. [l^. To investigate such finite size effects we have also 
performed simulations with nipsL < 3. 

We compute the quark propagators using point sources which we (Jacobi) smear to improve 
overlap with the ground state. For the three-point correlation functions we apply standard sequen- 
tial source techniques. The distance between source and sink is about 1 fm. Throughout this paper 
we will ignore contributions coming from disconnected terms. While these anyhow cancel in the 
iso-vector cases, results for the iso-scalar case maybe affected by an uncontrolled systematic error. 

To set the scale we use the Sommer parameter ro/a which we extrapolated to the chiral limit at 
each beta. While on the lattice this quantity can be determined with small statistical errors, there is 
no experimental determination. We therefore computed the dimensionless quantity {amN){ro/a) on 
the lattice and use the experimentally well known mass of the nucleon niN to obtain = 0.467 fm. 
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Figure 1 : The left panel shows the simulation points in the mpg vs. m^^L plane. In the right panel 
dashed lines show the lattice spacing and box sizes for which simulations have been performed. In 
both figures the continuous lines show where mpsL = 3. 



Most of the quantities considered in this paper need to be renormalised. The renormalisation 
constants have been determined using the RI'— MOM scheme [^], except for the vector current 
renormalisation constant Zy. Here we applied the condition that the nucleon's local vector current 
at zero momentum must be I. If necessary, the results are converted into MS scheme using the 4- 
and 2-3-loop expressions of the j8 function and corresponding anomalous dimension y, respectively. 



3. Lowest moments of PDFs 

Let us first consider the lowest moment of the polarized nucleon PDF (I)a^ (also known as 
axial coupling constant gA)- This quantity is determined from the renormalised axial vector current 
= Za (1 + <^fnq)A^, where arriq = (1/k" — 1/ Kc^^)/2. Za is known non-perturbatively for 
bA we use a tadpole improved one-loop perturbation theory result. 

We have fitted our lattice results to an expression from ChEFT based on the SSE formalism. 
Using this formalism both the quark mass dependence and the finite volume dependence [||] 
have been calculated. Since our results for different lattice spacings do not exhibit clear discretiza- 
tion effects we combine all our results where mps < 450 MeV. The fit range has been chosen such 
that stable fits are obtained. Our data is not sufficiently precise to determine all parameters. We 
therefore fix a few parameters to their phenomenological value and keep only g/^ in the chiral limit, 
the leading AA-coupling gi and the SSE coupling term Bg{X) as free fit parameters. The resulting 
fit and the lattice data are shown in Fig. |2^. 

In our fit we only included results for the largest lattice at a given set of bare parameters. For 
some data sets we have results for different volumes. We thus can compute the relative shift 

both from the fit as well as from the lattice data taking the results on the largest lattice as approxi- 
mation of gK{°°). In Fig. ^ we compare the relative shift for different values of the quark mass with 
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(a) (b) 
Figure 2: The left panel shows gA as a function of m^^. The open and filled diamonds show the 
lattice results before and after correction of finite size effects, respectively. The star indicates the 
experimental result. The line shows a fit to the data as described in the text. The right panel shows 
the relative finite size effects determined on the lattice (symbols) and obtained from a fit to an 
expression from ChEFT. 



our lattice results at mps ~ 270 Me V. The shift predicted from ChEFT only slightly underestimates 
the relative shift computed on the lattice. 

Also after correcting for finite size effects we observe a significant difference to the experi- 
mental value. It is interesting to notice that a much better agreement with the experimental value is 
observed for the ratio Ea/Ips (see Fig. In this ratio the renormalization constant Za drops out. 

In Fig. |3^ we show our results for the nucleon tensor charge (1)5^ = ,?t- We observe only a 
very mild quark mass dependence and the data reveals no systematic discretization effects. This 
quantity is not well known experimentally. Our values are larger than the phenomenological results 
presented in [^. 
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Figure 3: The left panel shows ^a/Zps as a function of mpg. The right panel shows our results for 
g^^ at a scale =2GeV. 
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4. n = 2 moments of PDFs 



The lowest moment of the unpolarized PDF {x)q = V2 corresponds to the momentum fraction 
carried by the quarks in the nucleon. Lattice results from different collaborations tend to be signifi- 
cantly larger than the phenomenological value. Fig. ^ and 4b show our most recent results for the 
iso-vector and iso-scalar channel. In the latter case disconnected contributions have been ignored. 

Also shown are the results from a fit to results utilizing methods of covariant Baryon Chi- 
ral Perturbation Theory (BChPT) [Q]. Fits have been performed with most parameters fixed to 
phenomenologically known values. The iso-vector (iso-scalar) channel data is fitted with only 2 
free parameters: V2 in the chiral limit and the coupling cs (cg). Near the physical light quark 
masses, BChPT predicts V2 to become larger when the quark mass becomes heavier In our data for 
mps < 250 MeV we do not see any indication for a bending down when approaching the physical 
pion mass. It thus does not seem that a lack of results at sufficiently small quark masses could 
explain the large discrepancy between the phenomenological value and the lattice results. There 
are some indications that part of the discrepancy can be explained by excited state contamination 
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Figure 4: The left and right panel show results for the second moment of the iso-vector and iso- 
scalar unpolarized PDFs, respectively, as a function of nip^. The solid lines show the fits to an 
expression from ChEFT. 



In Fig. 5a the results for the second moment of the polarized PDF {x)^q = a\ is shown. Dis- 
cretization effects again seem to be absent in data. From a comparison of the results for different 
volumes it seems that also finite size effects are small. Results from Heavy Baryon Chiral Pertur- 
bation Theory (HBChPT) [S] lead to the following expression: 
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(4.1) 



In Fig. ^ we plot this expression using pL = hzn and C chosen such that it matches the phenomeno- 
logical value. The bending down which we observe in our data for mps < 0.5 MeV is much less 
than one would expect from this HBChPT result. 
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Figure 5: The left panel shows the results for the second of the polarized PDFs as a function of 
. In the right panel the results for the Dirac form factor radius ri are plotted. The dashed lines 
show results from ChEFT as described in the text. 



5. Electromagnetic form factors 



To compute the electromagnetic form factors one makes use of the standard decomposition of 
the nucleon electromagnetic matrix elements {p' ,s'\V^\p,s) =u Y^Fi{Q^) + ^^^F2{Q^) u, (in 
Euclidian space) where we use the local vector current . p {s) and p' {s') denote initial and final 
momenta (spins), q = p' — p the momentum transfer (with = —q^). To calculate form factor 
radii and the anomalous magnetic we have to parametrize the lattice results. Here we use the ansatz 



Fi{Q' 



1 + 



(5.1) 



with p = 2 and p = 3 for the Dirac and Pauli form factors F\ and F2, respectively. Our data is not 
sufficiently precise to favour a particular parametrization (see [ |To| ] for another parametrization). 



From fits to Eq. (5.1 ) we determine the form factor radii r\ and r2 as well as the anomalous 
magnetic moment K. The quark mass dependence of these quantities has been calculated using 



the SSE formalism [11]. For r\ the parameters are known and we therefore restrict ourselves to a 



comparison of the SSE result and the lattice data (see Fig. 5b). While for mps > 300MeV the lattice 
results are significantly smaller than the phenomenological value, towards smaller quark masses we 
observe an increase of the radius. This is consistent with predictions from ChEFT For rj and K we 
find a similar behaviour. Since there are no phenomenological estimates for all parameters of the 
SSE expressions we perform a combined fit. The results are plotted in 6a and BR 



6. Summary and outlook 

We have presented an update of QCDSF results on the lowest moments of unpolarized, po- 
larized and tensor PDFs as well as the electromagnetic form factors. Some of our results at light 
quark masses with mps < 300MeV confirm the expectations from ChEFT that light quark mass 
effects are significant. However, this possibly does not explain all of the observed discrepancies 
from phenomenological values. 
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Figure 6: The left panel and right panel shows the results for the Pauli radius r2 and the anomalous 
magnetic moment k. The solid lines show fits to an expression from ChEFT. 
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